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^_' ■ Abstract. Two random processes X and Y on a metric space are said to be e-shy coupled 

(<— ^ I if there is positive probability of them staying at least a positive distance e apart from each 

^S| ■ other forever. Interest in the literature centres on nonexistence results subject to topological 

5^ I and geometric conditions; motivation arises from the desire to gain a better understanding 

^ • of probabilistic coupling. Previous non-existence results for co-adapted shy coupling of 

reflected Brownian motion required convexity conditions; we remove these conditions by 

showing the non-existence of shy co-adapted couplings of reflecting Brownian motion in any 

bounded CAT(O) domain with boundary satisfying uniform exterior sphere and interior cone 

conditions, for example, simply-connected bounded planar domains with C^ boundary. 

The proof uses a Cameron-Martin-Girsanov argument, together with a continuity prop- 

Pm ■ erty of the Skorokhod transformation and properties of the intrinsic metric of the domain. 

Ph I To this end, a generalization of Gauss' Lemma is established that shows differentiability 

r^ ■ of the intrinsic distance function for closures of CAT(O) domains with boundaries satisfy- 

"t^ I ing uniform exterior sphere and interior cone conditions. By this means, the shy coupling 

question is converted into a Lion and Man pursuit-evasion problem. 



cn 



1. Introduction 



en 
> 

,__! ! 1.1. Results and motivation. Benjainini, Burdzy, and Chen (2007) introduced the no- 

cn I tion of shy coupling: a couphng of Brownian motions X and Y (more generally, of two 

t~^ ■ random processes X and F on a metric space) is said to be shy if there is an e > such 

O , that P [dist(X(t), Y{t)) > e for all t] > 0. For example consider Brownian motion X on the 

circle: if Y is produced from X by a non-trivial rotation then X and Y exhibit a shy coup- 
ling, since dist{X,Y) is then constant. Interest in the existence or non-existence of such 
rj, - couplings arises from the study of couplings of reflected Brownian motions, which occur 

%j ■ in various contexts. Benjamini et al. (2007) discussed existence and non-existence of shy 

5^ I couplings for Brownian motions on graphs and for reflected Brownian motions in domains 
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(connected open subsets of Euclidean space) satisfying suitable boundary regularity condi- 
tions. They restricted attention to Markovian couplings and we will do essentially the same, 
by restricting attention to co-adapted couplings. (This is only slightly more general, but 
is more convenient for expression in terms of stochastic calculus.) In particular the results 
in Benjamini et al. (2007) showed that no shy co-adapted couplings can exist for reflected 
Brownian motion in convex bounded planar domains with C^ boundary satisfying a strict 
convexity condition (namely, that the boundary contains no line segments). Their argument 
used a large deviations argument bearing some resemblance to methods from differential 
game theory. Kendall (2009) showed that neither differentiability nor strict convexity is 
required for the planar result, and also generalized the result to convex bounded domains in 
higher dimensions whose boundaries need no longer be smooth but still satisfy the regularity 
condition requiring triviality of all line segments contained in the boundary. These more 
recent results are based on direct proofs using ideas from stochastic control. 

The work described below both generalizes the above results and also shows that absence 
of shyness is not confined to the case of convexity. We consider a bounded domain with 
boundary satisfying uniform exterior sphere and interior cone conditions and that satisfies 
a CAT(O) condition (see Definition 4) when furnished with the intrinsic metric, and we 
show that such domains cannot support shy co-adapted couplings of reflected Brownian 
motions. We do this by establishing a rather direct connection between (the non-existence 
of) Brownian shy co-adapted couplings and deterministic pursuit-evasion problems. As part 
of this process, we generalize Gauss' lemma (on the differentiability of the distance function) 
to the case of closures of CAT(O) domains furnished with the intrinsic metric and satisfying 
uniform exterior sphere and interior cone conditions. It may not be evident to the reader 
exactly how the stochastic and undirected notion of Brownian motion can be connected to 
the deterministic and intentional notion of a pursuit-evasion problem, and it was not initially 
evident to us (though, in retrospect, this is latent in Benjamini et al., 2007), but nonetheless 
the connection is both immediate and useful. 

The pursuit-evasion problem in question is a well-known problem concerning a Lion chas- 
ing a Man in a disk, both travelling at unit speed: R. Rado's celebrated "Lion and Man" 
problem. Our shy coupling problem leads us to consider the generalization in which the Lion 
chases the Man in a bounded domain which is CAT(O) in its intrinsic metric. Isaacs (1965) 
is the classic reference for pursuit-evasion problems; Nahin (2007) provides an accessible 
exposition of the special case of the Lion and Man problem in the unit disk. Littlewood 
(1986, pp. 114-117 in Bollobas' extended edition) provides a brief description of the Lion 
and Man problem with an indication of its history, including a presentation of Besicov- 
itch's celebrated proof that in the disc the Man can evade the Lion indefinitely, even though 
the distance between Lion and Man may tend to zero. A generalization of discrete-time 
pursuit-evasion to bounded CAT(O) domains is dealt with in Alexander, Bishop, and Christ 
(2006); we summarize concepts from metric geometry and develop results required for the 
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continuous-time variant in Section 2, and it is here that we generahze the Gauss lemma to 
the case of closures of CAT(O) domains with sufficient boundary regularity (Proposition 14). 

In particular, Section 2 rigorously develops the geometric results required to reason with 
these concepts in the context of the intrinsic metric for the domain D (determined by lengths 
of paths restricted to lie within D). On a first reading one should feel free to note only the 
general ideas of Section 2, and then to pass quickly on to the probabilistic arguments in the 
remaining sections of the paper. 

In Section 3, we describe how continuous-time pursuit-evasion problems can be solved 
in CAT(O) domains. We obtain an upper bound for the time of e-capture, expressed in 
terms of domain geometry. Simultaneously with and independently of our research project, 
Chanyoung Jun developed in his Ph.D. thesis (Jun, 2011) a theory of continuous pursuit in 
CAT(k) spaces that overlaps somewhat with our results. 

Pursuit-evasion games involve control of the velocity of the pursuer so as to bring it 
arbitrarily close to the evader, regardless of what strategy may be adopted by the evader. 
In order to show non-existence of Brownian shy couplings, we investigate the possibility of 
bringing the Brownian pursuer (the Brownian Lion) arbitrarily close to the Brownian evader 
(the Brownian Man), regardless of how the Brownian motion of the Brownian Man is coupled 
to that of the Brownian Lion. The connection between coupling and deterministic Lion and 
Man problems is described in Section 4: a suitable pursuit strategy generates a vectorfield 
X on the configuration manifold generated by the locations of Brownian Lion and Man. 
(More pedantically, it generates a section of the puUback of the tangent bundle of D to the 
configuration space of the pursuer and evader before capture.) If this pursuit strategy can 
be guaranteed to bring the Lion within e/2 of Man by a bounded time t^. in the deterministic 
problem, then a Cameron-Martin-Girsanov argument together with a continuity property 
for the Skorokhod transformation shows that the Brownian Lion has a positive probability 
of getting within distance e of the Brownian Man, whatever coupling strategy might be 
adopted by the Brownian Man. 

The paper concludes with Section 5, which discusses possible extensions of these results, 
further questions, and conjectures. 

We now state the main results of this paper, using terms defined in Section 2. Here and 
elsewhere in the paper, we consider only domains in Euclidean space of dimensions 2 or 
higher. 

Theorem 1. Suppose that D is a hounded domain with boundary satisfying uniform exterior 
sphere and interior cone conditions, and which is CAT{{]) in its intrinsic metric. There can 
he no shy co-adapted coupling for reflected Brownian motion in D. 

Examples of CAT(O) domains include convex domains and domains that are the unions 
of a pair of convex domains. See, for instance, Bridson and Haefiiger (1999) and Alexander 
et al. (2006), where more general examples are also provided; in particular, a large range of 
examples follows from iterated application of the result that if two CAT(O) domains have a 
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geodesically convex intersection then their union is CAT(O). The exterior sphere and interior 
cone conditions in the theorem are required in order to apply the results of Saisho (1987) to 
generate reflected diffusions using the Skorokhod transformation. 

The three dimensional domain in Figure 1 is CAT(O). There are two different ways to see 
this. First, it is easy to see that for every point on the boundary of the domain, at most one 
of the principal curvatures is negative. An alternative way to see that the domain is CAT(O) 
is to observe that a single dumbbell (the union of two spheres and the connecting tube) is a 
CAT(O) domain. The whole set is the union of five dumbbells. The non-empty intersections 
of the dumbbells are balls. 




Figure 1. A CAT(O) example which is the union of five dumbbells. 



Remarkably, all bounded simply-connected planar domains are CAT(O) in their intrinsic 
metrics. Thus, in the planar case, there is an immediate consequence of Theorem 1 which is 
a strikingly powerful result depending principally on topological conditions: 

Theorem 2. Suppose that D is a simply- connected hounded planar domain with boundary 
satisfying uniform exterior sphere and interior cone conditions. There can he no shy co- 
adapted coupling for reflected Brownian motion in D. 

1.2. Some basic tools for probabilistic coupling. All probabilistic couplings considered 
here are co-adapted couplings, which are defined for general Markov processes in Kendall 
(2009). In essence, a co-adapted coupling of two Markov processes is a construction of 
the two Markov processes on the same probability space, which are adapted to the same 
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filtration such that each process possesses the prescribed transition functions with respect to 
the common filtration. 

In this paper, it suffices to work with co-adapted couplings of d- dimensional Brownian 
motions: B and B are said to be co-adaptively coupled Brownian motions if they are defined 
on the same probability space and adapted to the same filtration {J-'j : t > 0} and if, in 
addition, both satisfy an independent increments property taken with respect to the common 
filtration: 

Bt+s — Bt is independent of J-t for alH, s > , 
Bt+s — Bt is independent of J-t for alH, s > . 

Note that -Bf+s — Bt and -Bj+s — Bt need not be independent of each other. Kendall (2009, 
Lemma 6) shows that one may represent such a coupling using stochastic calculus, possibly 
at the cost of augmenting the filtration by adding a further independent Brownian motion 
C: there exist {d x (i)-matrix-valued predictable random processes J and K such that 

B = /"j^dfi+ [K^dC- 



moreover, one may choose J'''J + ]K''^K to be equal to the {d x d) identity matrix at all times. 
A pair of processes X and X is said to form a co-adapted coupling if they can be defined 
by strong solutions of stochastic differential equations driven hj B, B respectively. In the 
paper, we will employ the stochastic differential equation obtained from the Skorokhod 
transformation for reflected Brownian motion in a domain D of suitable boundary regularity, 
such as under uniform exterior sphere and uniform interior cone conditions, as discussed in 
Section 2. For r > 0, set N'x,r = {z^ € M*^ : |z/| = l,B{x + rUjr) (1 D = 0}. The vectors z/ 
can be be viewed as "exterior normal unit vectors at x G dD" ; note that there may be more 
than one such vector at a particular point x G dD. The set A^,,- is decreasing in r, and 
the uniform exterior sphere condition asserts that r can be chosen so that, for all x G dD, 
■f^x,r 7^ 0, with N'x,r = ■^x,s for < s < r. Under uniform exterior sphere and uniform 
interior cone conditions, Saisho (1987) has shown that, given a driving Brownian motion B, 
there exists a unique solution pair (X, L^) satisfying 

dX = AB-vxdL^, 

L is non-decreasing and increases only when X G dD , 

VX e A/'x,r • 

Thus L^ may be viewed as the local time of the reflected Brownian motion X on the 
boundary dD. 

In this paper, all vectors are assumed to be column vectors unless specified otherwise. 
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2. CAT(O) GEOMETRY AND THE DETERMINISTIC PURSUIT- EVASION PROBLEM 

Recall that the intrinsic metric for a domain D is generated by the infimum of Euclidean 
lengths len(7) of smooth connecting paths 7 lying wholly within the domain. (The definition 
is typically formulated in the context of general metric spaces and regularizable paths.) 

Definition 3. The intrinsic distance between two points x and y in a domain D is given by 

(1) distintr{x,y) = inf {len(7) : 'j is a smooth path connecting x and y in D} . 

For a domain D, a standard compactness argument shows that paths attaining the infimum 
of (1) will always exist in the closure of the domain: these are called intrinsic geodesies. 

As described in Bridson and Haefiiger (1999, §11.1, Definition 1.1) (see also Burago, Bur- 
ago, and Ivanov, 2001), one can define simple curvature conditions for metric spaces such as 
{D, distintr), based on the behaviour of geodesic triangles. We first give the case of comparison 
with fiat Euclidean space (which has zero curvature). 

Definition 4. We say that {D, distintr) is a CAT{0) domain if the following triangle com- 
parison holds: Suppose that Ta^, Ta^c, (md F^c are unit-speed intrinsic geodesies for D, 
connecting points a to b, a to c, and b to c respectively. Then, for all such geodesic triangles, 

dist intr(Xa,b{^),'ra,c(t)) < r{s,t) , 

where r{s,t) is the distance between points at distance s, respectively, t, from a along the 
side ab, respectively ac, of an ordinary Euclidean triangle abc that has the same side-lengths. 





Figure 2. Illustration of the CAT(O) condition. 

Consequently, chords of triangles in {D, distintr) are shorter than comparable chords of the 
comparable Euclidean triangles, as illustrated in Figure 2. 

Bridson and Haefiiger (1999, §11.1, Definition 1.1) actually introduces the more general 
notion of a CAT(/i;) domain (see also Alexander, Bishop, and Christ, 2010, Appendix A). 
Here we describe the case when comparisons are drawn with triangles on a sphere of radius 
1/a/k, for K > (hence the sphere has curvature k). It is necessary here to restrict attention 
to suitably small triangles, as measured by perimeter. 
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Definition 5. We say that D is a CAT{k,) domain for k, > if any two distinct points with 
distance less than tt/a/k are joined by a geodesic and the distance between any two points 
of any geodesic triangle Apqr of perimeter less than 27r/i/K is no greater than the distance 
between the corresponding points of the model triangle Apqr with the same sidelengths in the 
2-dimensional Euclidean sphere of radius 1/\/k. 

Remark. Gromov introduced the acronym CAT, standing for Cartan, Aleksandrov, Toponogov. 
In this paper, we will mostly consider spaces CAT{k,) with n = 0. We include some results 
concerning the CAT{k) case with k, > because they will be used in the forthcoming paper 
Bramson, Burdzy, and Kendall (2011). 

Remark. ^45 noted in Bridson and Haefliger (1999, Proposition II. 3.1) (see also Burago 
et al., 2001, %4-3), in CAT{k) spaces the notion of angle is well-defined for (locally) minimal 
geodesies. 

Consequently, geodesies in a CAT(O) space diverge at most as fast as corresponding 
geodesies in Euclidean space. Note that CAT(O) is a global condition, applying to all 
possible geodesic triangles. In particular it can be shown that CAT(O) spaces are always 
simply-connected and indeed contractible (Bridson and Haefliger, 1999, Proposition II. 1.4, 
or Alexander et al., 2010, Appendix A). 

Remarkably, bounded planar domains are CAT(O) if they are simply-connected; see Bishop 
(2008) for a careful proof. Readers may convince themselves of this at an intuitive level by 
drawing pictures (as exemplified in Figure 3); as is the case with other foundational results 
in metric spaces, the rigorous proof requires delicate reasoning. 





Figure 3. Illustration of the CAT(O) property for a bounded simply- 
connected planar domain. The effect of the boundary is to make the triangle 
"skinnier" than its Euclidean counterpart, thus establishing the CAT(O) com- 
parison property. 

We now introduce two complementary notions of boundary regularity following Saisho 
(1987). An exterior sphere condition (also called weak convexity) requires that every bound- 
ary point is touched by at least one external sphere. Here and in the following, let B{y^ s) 
denote the open Euclidean ball of radius s centered on y. 
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Definition 6 (Uniform exterior sphere condition, from Saisho, 1987, §1, Condition (A)). 
A domain D is said to satisfy a uniform exterior sphere condition, based on radius r if, for 
every x G dD, the set of '^exterior normals" Mx,r = {/^ G M'^ : |z/| = 1, B{x + rv^ r) fl Z^ = 0} 
is non-empty, with N'x,r = A4,s for < s < r. 

Thus a uniform exterior sphere condition allows one to move a fixed ball all the way 
around the outside of the domain boundary. In particular, D can have no "inward-pointing 
corners". Here is a simple observation which will be useful later and corresponds to the 
intuition about being able to move a fixed ball about D; such D may be represented as 
intersections of complements of balls, in a manner entirely analogous to the representation 
of a convex set as the intersection of half-planes (so justifying the alternative term "weak 
convexity" ) . 

Lemma 7. Suppose that the domain D satisfies a uniform exterior sphere condition based 
on radius r. Then 



D = f]{B{z,ry : i3(^,r)nD = 0} 



Proof. Let the Minkowski sum A (B B oi two Euclidean sets A and B he A (B B = {x + y : 
X E A,y E B}. Certainly F = f]{B{z,rY : B{z,r) fl D = 0} is closed, since B{z,r) is an 
open ball. Moreover D (1 F; hence D 'O F. Furthermore F C D © B{o,r), where o is the 
origin of the ambient Euclidean space. 

Following Saisho (1987, Remark 1.3), because of the uniform exterior sphere condition, we 
can define a projection x i— ?■ x from D(BB{o, r) onto D using the Euclidean metric. Consider 
X E D (B B{o,r). Then the projection x E D is defined; moreover, if x G dD and x ^ D, 
then 

X — X 



A4. 



\x — x\ 



is a unit vector whose offset produces a tangent sphere of radius r at x (using the argument 
of Saisho, 1987). But this implies that if x G (D®B{o,r)) \ D then 

„ /_ X - X 
X G o X — r— -,r 

\ |x — x| 

and so X ^ F. Accordingly D = F as required. D 

On the other hand, a uniform interior cone condition requires that any boundary point 
supports a bounded cone truncated to the boundary of a ball, and moreover that the cone 
may be translated locally within the domain. 

Definition 8 (Uniform interior cone condition, from Saisho, 1987, §1, Condition {B')). A 
domain D is said to satisfy a uniform interior cone condition, based on radius S > and 
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angle a G (0,7r/2], if, for every x G dD, there is at least one unit vector m such that the 
cone C(m) = {z : {z, m) > 1^1 cos a} satisfies 

{y + C{m))nB{x,6) C D for all y e D n B{x,6) . 

We say that the cone y + C(m) is based on y and angle a G (0, vr/2]. 

Thus a uniform interior cone condition implies that the "outward-pointing corners" must 
not be too sharp. Note that Saisho actually uses a slightly weaker condition with less intuitive 
content (Saisho, 1987, Condition (B)); we do not consider this weaker notion further in what 
follows. 

In fact, the property of a domain satisfying a uniform interior cone condition is equivalent 
to it being a Lipschitz domain. 

Definition 9 (Lipschitz domain). Recall that a function f : M'^"^ -^ M is Lipschitz, with 
constant A < oo, if \f{x) — f{y)\ < A|a; — y\ for all x,y E W^~^ . A domain D is said to he 
Lipschitz, with constant A, if there exists 5 > such that, for every x G dD, there exists 
an orthonormal basis Ci, 62, • • • , e^ and a Lipschitz function f : W^~^ — )■ M, with constant X, 
such that 

Bix,S)nD = {yeBix,6):fiyi,...,yd^i)<yd}, 
where we write yi = (y^ei) , . . . , yd= {y, e^). 

The equivalence of Definitions 8 and 9 depends on the fact that the cone axis vector in 
Definition 8 is chosen to be the same for all y G B{x,6), and so can be used as e^ in the 
orthonormal basis for B{x, 6) required in Definition 9. The constants A and a in Definitions 
8 and 9 are related by A = cot a, while the two 5's of Definitions 8 and 9 may be taken to 
be equal. Note too that if the uniform interior cone / Lipschitz domain property holds for a 
given 6 > 0, then evidently it also holds for all smaller 6. 

If a domain satisfies a uniform interior cone condition, then the intrinsic metric and Euc- 
lidean metric properties are closely related. 

Lemma 10. A domain D that is bounded in Euclidean metric and satisfies a uniform interior 
cone condition must have finite intrinsic diameter. 

Proof. Certainly distintr(3^, y) is a continuous function of {x,y) in the open set D x D and 
takes only finite values there. Note that the domain D is path-connected, being an open 
connected subset of Euclidean space. 

Suppose that D satisfies a uniform interior cone condition based on radius 6 > and angle 
a G (0, 7r/2]. If m is a unit vector for the interior cone condition at a; G dD then geometrical 
arguments show that x + |5m is at least ^Ssina from the exterior D'^. Choosing 6' with 
< 6' < l^sina, it follows that any such x + |5m belongs to 

DeB{o,6') = {D'®B{o,6')Y = {x e D : B{x,6') C D} , 
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which itself is closed. Inheriting boundedness from D, it is therefore compact in the Euclidean 
topology, and hence also in the topology derived from the intrinsic metric, since the two 
metrics are locally equal away from the boundary of D. Hence {distintr(2^,|/) : x,y E D Q 
S(o, 5')} attains a maximum value, which is therefore finite. However, for any x' ,y' G -D, 
we have 

(2) distintr(a;', y') < S + sup{distintr(a;, y) : x,y e DQ i3(o, S')} , 

because we have used the uniform interior cone condition to ensure that from each point 
on the boundary there is a straight-line segment of length |(5 to {distintT{x,y) : x,y E 
D Q i3(o, (5')}. Hence the intrinsic diameter must be bounded by the right-hand side of 

(2). D 

The full force of the uniform interior cone condition is not required for the above result; the 
proof does not require coordination of the directions of interior cones at different base-points. 
The full force of the uniform interior cone condition assures us that any path of finite length 
leading in D to a point x on the boundary of D can be deformed continuously in D into one 
which in its final phase is the segment on which the interior cone at x is based. Moreover, 
the lengths of the curves throughout this deformation can be constrained to be arbitrarily 
close to the length of the original path. This allows us to view D as a topological manifold 
with boundary, which is continuously embedded in the ambient Euclidean space. More than 
this, it shows that the completion D oi D under the intrinsic metric can be identified with 
the Euclidean closure D and moreover that the intrinsic metric and the Euclidean metric 
actually endow D with the same topology. Finally, Bridson and Haefiiger (1999, Corollary 
II. 3. 11) show that the closure D, viewed as the completion D oi D in intrinsic metric, inherits 
CAT(O) structure from D. 

2.1. Regularity for geodesies. We wish to consider pursuit-evasion in a bounded CAT(O) 
domain. Lion and Man both move with unit speed, with the Lion seeking to draw closer 
to the Man by using a "greedy" pursuit strategy (which is not necessarily optimal). This 
Lion strategy can be phrased in terms of an M*^- valued field x of unit vectors defined on the 
configuration space {D x D)\ {{x, x) : x E D}, such that xi^^ v) is the initial velocity of the 
unit-speed geodesic moving from x to y. (This is the vectorfield described pedantically in 
Section 1 as a section of the puUback of the tangent bundle of D to the configuration space 
of the pursuer and evader before capture.) 

We first show that the combination of uniform exterior sphere and uniform interior cone / 
Lipschitz conditions implies that, working locally, every boundary point of the intersection 
of the domain D with a suitable 2-plane will support an exterior sphere, albeit with smaller 
radius. 

Lemma 11. Suppose that D is a domain satisfying a uniform exterior sphere condition 
based on radius r > 0, and a uniform interior cone condition based on radius 5 > and 
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Figure 4. Illustration of interior cone C^ and exterior ball B{y,r) at 2 G -D. 

angle a G (0,7r/2]. Suppose that z G dD and e^ is the d*^ vector in the orthonormal basis 
corresponding to z as in Definition 9. Let P he a 2-plane intersecting D and containing 
z and z + e^. Then there exists w G P, with \w — z\ = dist(tf;,-D fl P) = rsina and 
B{w,rsma) n (D n P) = 0. 

Since the interior cone condition is uniform, the lemma shows that every point in the 
boundary of D fl P near z must support an exterior sphere of radius r sin a. 

Proof. Suppose that z G dD with e^ defined as above. Let P be a 2-plane containing z and 
z + Cd- Since z G dD, there is an exterior sphere touching z, defined by a ball B{y,r) C D"^ 
with z G B{y, r). By Definition 9, the cone 

Cz = {w : {w — z^Cd) > \w — z\cosa} 

lies locally in D, in the sense that C^ H B{z,5) C D (see Figure 4). If | + /3 is the angle 
between Cd and y — z, then two-dimensional geometry (Figure 5) shows that 

min {|i/ — (76^ + 2;)| : 7 G M} = rcos/3. 

But f3 > a ii C^n B{z, 6) C D and B{y, r) C D'^] moreover, the line {76^ + 2 : 7 G M} must 
lie in P. Hence the distance from y to P is at most r cos P < r cos a. Consequently the 
radius of the disk B{y,r) n P is at least rsina; since z G d{D fl P) and B{y,r) fl P is an 
exterior sphere to z in P, the lemma follows. D 

We can now establish some important technical consequences of the uniform exterior 
sphere and interior cone conditions together with the CAT(O) condition; namely, that the 
Euclidean and intrinsic distances are locally comparable, and that the vectorfield x is con- 
tinuous with reference to the common topology of the Euclidean metric and the intrinsic 
metric, and hence is uniformly continuous over regions for which the two arguments are 
well-separated. This is spelled out in the following proposition. In fact we state and prove a 
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'B(y,r) 



Figure 5. Two-dimensional section of Figure 4 illustrating the underlying 
two-dimensional geometry. 

generalization of the result for CAT{k,) domains with k > 0, so that we can apply it in the 
forthcoming paper Bramson et al. (2011). 

Proposition 12. Suppose that D is a CAT{k,) domain with k > 0, bounded in the Euclidean 
metric and satisfying a uniform exterior sphere condition based on radius r > 0, and a 
uniform interior cone condition based on radius 6 > and angle a G (0,7r/2]. We can and 
will assume without loss of generality that A = cot a > 1 . 

(1) Suppose a, b & D are close in the Euclidean metric, in the sense that 

(3) |a -6| < min{(5/(4A),2rsina}. 

Then 

distintr{a,b)' 



(4) 



2rsina sin 



< 



b\ < distintria, b) 



2r sin a 
0. Intrinsic geodesies for D (necessarily minimal, by the CAT{0) 



(2) Suppose that k 

condition) are continuously differentiable and their direction fields satisfy a Lipschitz 
property with constant A= ^^ \^ ^ that therefore holds uniformly for all minimal in- 



trinsic geodesies in D and hence in D (since CAT{^) geodesies depend continuously 
on their endpoints). For k > 0, the same conclusion holds for minimal intrinsic 
geodesies with endpoints in D which are separated by intrinsic distance strictly less 
than ir/y/Ji. 
(3) For X, y in D with distintr{x, y) < t^ j \fi^ (as usual, tt/v^ = oo^, let x(x, y) be the unit 
vector at x pointing along the unique intrinsic geodesic ^(^'S') from x to y. Then x(x, y) 
depends continuously on {x,y) in A = {{x,y) E D x D : < distintr{x,y) < 7r/y^} 
and hence is uniformly continuous over compact subregions of A. 
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Proof of Proposition 12, part (1). Definitions 8 and 9 are equivalent, so the domain D is 
Lipschitz with constant A = cot a. For each ball of radius 5, we may therefore construct a 
coordinate system ei, . . . , e^ and a Lipschitz function / to implement the Lipschitz property 
oiD. _ 

Consider a, b ^ D with |a — 6| < min{(5/(4A), 2r sina}. If the line segment 5* between 
a and b does not intersect dD, then it must form the (unique, minimal) intrinsic geodesic 
between a and b, and (4) follows immediately. If 5* does not intersect mt{D'^), then we can 
cover the intersection 5* fl dD with a single B{z, 5) (for z G dD) and use the unit vector 
Cd corresponding to the ball (equivalently, the unit vector defining the cone for the ball) to 
perturb 5 to a regularizable path in D (save for the endpoints) with length arbitrarily close 
to that of S. Hence S is the intrinsic geodesic between a and b, and therefore (4) follows 
immediately. So we can confine our attention to the case when a ^ b and S intersects 
int(L'^). 

Applying Definition 9 to B{a, 6), there is a Lipschitz function / : M'^"^ — )■ M, with Lipschitz 
constant A = cot a, and an orthonormal basis ei, . . . , e^, such that 

(5) B{a,6)nD = {yeB{a,6):f{yi,...,yd-i)<y4 
where yi = {y,ei), . . . , y^ = {y, e^). Consider 

C{a) = {yeR": (|yi - ai\^ + . . . + \yd-i - a^-ilY^ < 5/{A\) , \yd - a^] < 6/2} 

and note that, since A > 1, it is a consequence of (5) that C(a) C B{a,6). Moreover / has 
Lipschitz constant A, so we can control the behaviour of that part of the boundary of D 
lying within C{a): 

(6) {yeR": (|yi - a,\^ + . . . + ly^-i - a^^il^' < 6/{4X), f{y,, . . . , y^-i) = yd] 

= {ye C{a) : /(yi, . . . , |/,_i) = yd} = dD n C{a) . 

Applying Lemma 11 to the 2-plane 

P = a + linear span {b — a, e^} , 

every boundary point of Z^ fl P fl C{a) supports an exterior disk of radius r sin a. (Note that 
the Lipschitz representation implies that (dD) fl P fl C{a) = d{D fl P) fl C{a).) We shall use 
these exterior disks to construct a short path between a and b. 

It follows from (3) and (6) that the two rays from a and b along the direction e^ must lie 
in P n D until they leave C{a): 

C(a)n{a + 7ed :7 > 0} C PnD, 

(7) C{a) n {6 + 76^ : 7 > 0} C P n P> . 

We set u, V to be the intersections of these rays with dC{a). For each rj G (0, 1), consider 
the point rja + {1 — ri)b and the open segment which is the intersection of the corresponding 



SHY COUPLINGS, CAT(O) SPACES, AND THE LION AND MAN 



14 





u 






V 


D 


i 










^d 






k^ 








M^^ ~~^-., 




Ja 




X '' ^ 




V / 


C(a) 





Figure 6. Illustration of proof of Proposition 12, part (1). The aspect ratio 
is not realistic — the height of C{a) must be at least twice its horizontal dia- 
meter. 



ray with C(a), namely 

C{a) n {rja + (1 - r/)6 + 76^ : 7 > 0} . 
It follows from (3) and (6) that a non-empty final sub-segment 

C{a) n {rja + (1 - r/)6 + 76^ : 7 > 7^} 

must lie m D r\ P. But then any exterior disk for 77a + (1 — 77)6 + 7^6^ has to avoid the rays 
defined in (7) as well as the above non-empty final sub-segments; it must not intersect the 
segments [a,u\ and [fc, w], and also may not intersect that portion of dC{a) which intersects 
rays {rja + (1 — vi)h + •jCd : 7 > 0} (see Figure 6). Consequently (since |a — &| < 2r sina) 
such an exterior disk must have center lying on the side of the line through a and b which 
is opposite to the side containing u and v, and must not intersect the complement of the 
segment S* in the line through a and b. 

The envelope of the boundaries of all such disks of radius r sin a in D fl P is formed by 
the complement of the segment S in the line through a and b together with the minor arc 
A of the circle of radius r sin a running through a and b. We can use A to generate a short 
path between a and 6 in D as follows. If A does not intersect one of the rays in (7) then A 
itself suffices; otherwise a still shorter path may be formed which lies wholly in D by making 
a short-cut using the relevant ray. In any case a small perturbation of A or the short-cut 
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version, using the vector e^, will provide a path in D from a to 6 of length less than the 
length of A plus an arbitrarily small increment. Calculation of the length of the minor arc 
A now leads to the desired bounds on distintr(fl, b) as given in (4). 

Proof of Proposition 12, part (2). Consider points a, h and c in D, lying in this order along 
an intrinsic geodesic F in D. We will need the geodesic F to be minimal. This is immediate 
in case k = 0; in the case /t > it follows if we require that the length of F is strictly less 
than vt/a/k. For some positive t < min{(5/(4A), 2r sina}, suppose that the intrinsic distances 
between a and b and between b and c are both equal to t. Since F is a minimal geodesic, the 
intrinsic distance between a and c must be 2t. Let pi = \a — b\, p2 = \b — c\ and p-^ = \a — c\ 
be the Euclidean distances between these three pairs of points and let vr — 6* be the interior 
angle at b in the Euclidean triangle abc. By the cosine formula, 

cos^ = -cos(.-^) = -P' + P'-P' = Pi- Pi- Pi, 

2piP2 2pip2 

The upper bound on t means we can apply (4) to the intrinsic and Euclidean distances 
between a, b and c. Hence pi <t, p2 <t and 

f 2t \ 

Pa > 2r sin a sin > 

\ 2r sin a J 

where the last step uses sin a > a — a^/6 if a > 0. Together with the cosine formula, these 
bounds for pi, p2 and ps yield 

cos^ > ^ — ^ ^^ 




2piP2 

\ 6 \ r sin a/ I I I 1/t 



> ^ — ^ ^ — = 2 1 



hence 



2t^ \ Q \r sin a 



9 I f t 

cos - > 1 




2 6 \ r sin a ^ 

Considering t < min{(5/(4A), 2r sina}, it follows by calculus that there exists a c(t) tending 
to zero with t such that 

8 9 < -^ U + c^ 

a/3 2r sin a \ 2rsina^ 

Suppose now that the intrinsic geodesic F has total length K. For any positive integer 
m > i^/min{5/(4A), 2r sina}, let Oq = a, Oi, . . . , Om-i, a^ = b he m + 1 points equally 
spaced along the geodesic, so that distintr(oj-i, ctj) = t for j = 1, . . . , m. Define gm '■ [0, K] — >■ 
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Mf^ to be the piecewise-linear curve interpolating gmiJK/m) = aj for j = 0, . . . ,m. By (8), 
all the angles between successive line-segments of the trajectory of Qm are bounded above by 

a/3 2r sin a \ 2rsina^ 

Define the directional unit vector field of the curve Qm by uJm{s) = gmi^)/\9m{^)\ for s where 
Qmis) is linear, and extend to all s using left-limits for s > and the right-limit for s = 0. 
Then, by the triangle inequality, 

/ N / M 4 t ( , . t \ flso — Si\ 

-^3 2r sm a \ zrsinaj \ t 



2rsina / t , i / / m 

sm I :^ — : 1 < l9m(s)| < 1, 



From (4), 

2r 

t \^2rsina 

hence we obtain the inequality 

(9) 

\9LiS2) - 9'miSl)\ < \uJmiS2) - W^(si)| + |a;^(s2) - g'm{s2)\ + Wmi^l) - Wm(Sl)| 

4 t f ,, t \/|s2-si| \ / 2rsina / t 

< -^ 1 + c(t) ^— + 1+21 sin 

y/^zrsma \ zrsinaj \ t J \ t \2rsma^ 

from which there follows a uniform bound on the absolute variation of the g'^ functions. Thus 

we can apply Helly's selection theorem to deduce that g'^ will converge along a subsequence, 

both pointwise and locally in L^, to a continuous limit h. It is immediate that gm converges 

uniformly to F, and F must be almost everywhere differentiable with limit h = V. Moreover, 

from (9) (and bearing in mind that c(t) — )■ with t), we may deduce that the derivative F' 

is Lipschitz with constant 

4 1 

1/3 2r sin a ' 
and indeed that F is continuously differentiable. 

Proof of Proposition 12, part (3). As noted above, the CAT(k) property of D implies 
that all geodesies between points x and y satisfying < distintr(a^) y) < tt/v^ are unique and 
minimal. Consider {x,y), {xn,yn) ^ {{v,z) & DxD : < distintr('y, 2) < tt/^/k} with x„ — )■ x 
and ?/„—>■?/ in the Euclidean metric; taking subsequences we may suppose that x{^n,yn) 
converges to a limit. Part (2) of the lemma establishes the uniform Lipschitz property of the 
direction fields of all minimal geodesies in D so, by the Arzela-Ascoli theorem, we can find a 
subsequence (a;„,^, ynj such that the geodesies from a:„,^ to y„^ must converge to a curve from 
X to y whose direction field is the limit of the direction fields of these minimal geodesies; 
hence its direction at x must be linifc xl^n^, 2/n.fe)- By minimality of the geodesies from Xr^. 
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to Hn^. and taking limits, the length of the limiting curve can be no greater than that of the 
unique minimal geodesic from x to y; therefore the limiting curve must also be a minimal 
geodesic from x to y. By the CAT(/t:) property, the two minimal geodesies from x to y must 
therefore be equal, and therefore it follows that lim^ xl^n^, Vuk) = xi^^ u)- K follows that any 
subsequence of {xn,yn) -^ {x,y) (convergence in Euclidean metric) must possess a further 
subsequence for which lim^ xi^n^yUnk) = xi^y v): and therefore lim„ xi^n, Un) = x{^^ v) must 
hold. This establishes continuity of x with reference to the Euclidean metric. D 

Remark. Part (1) of Proposition 12 may he used to show that x(a, h) is Hdlder(\) in its 
second argument b when a and b are well-separated. We omit this argument, as the result is 
not used in this paper. 

Remark. Setting p = \a — b\ and t = dist intr{(i,b). Inequality (4) can be rewritten as 

2rsina f t \ p 

sin < - < 1 . 



d \2rsintty t 

The following is a trivial but useful consequence of the above estimates: for some Ci > 0, 
depending on D, and all a,b E D, with \x — y\ < ci, 

(10) \a — b\ < distintr{o,,b) < 2\a — b\. 

Moreover, since sin0 > — |0^ ^/ > 0, 

1 t^ p 

< - < 1. 



;il) 1 < - < 1 + C2t' < 1 + C3p' 



6 4r^ sin^ a t 

The last inequality and (10) imply that for some C2,C3 < oo, depending on 6,r and a, and 
for p < min{(5, 2r sin a}, 

t 

P 

Proposition 12 makes it possible to quantify the extent to which short intrinsic geodesies 
may be approximated by Euclidean segments. 

Corollary 13. Suppose the assumptions on Proposition 12 hold, and that T is a unit-speed 
intrinsic geodesic with intrinsic length t < min{(5/(4A), 2rsina}. Then 

|r(t)-r(o)-r(o)t| < ^ ^ 



3 2r sin a 

Proof. Set p = |r(t) — r(0)| equal to the Euclidean distance between the two end-points of 
F; then p is bounded above by the intrinsic length t. Let 9 be the angle between r'(0) and 

r(t)-r(o). 

Proposition 12, (2) tells us that P' is Lipschitz with constant -j^ 2rs\na ' Hence 

(12) (rw.r'(o)) = i-i|rw-r(o)r > i-i(^^)^ 
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and this integrates to 



(r(t)-r(o),r'(o)) > (i-^^^,-^)t. 



e 



Consequently 



|2 



|r(t) - r(o) - r'(o)tr = |r(t) - r(o)r + |r'(o)tr - 2 (r(t) - r(o), r'(o)t) 

e \ , 16 t^ 



< p" + r - 2 1 - - ^ . „ t' < 



9 4r2 sin a / 9 Ar"^ sin a 

The resuh follows by taking square roots. D 

At this point we revert to considering CAT(O) spaces only, since generalization of the 
following proofs to the CAT(k) case would extend the exposition. We recall Gauss' lemma 
from Riemannian geometry, that the exponential map is a radial isometry. Cheeger and Ebin 
(2008, Chapter 1, Section 2) observe that, for smooth Riemannian manifolds, it is equivalent 
to the assertion that the Riemannian distance distintr(2;,|/) is continuously differentiable in 
X when x ^ y and y does not lie in the cut-locus of x, with the gradient being given by 
the tangent of the geodesic running from y to x. Proposition 12 and Corollary 13 can 
be used to prove the following Gauss lemma for CAT(O) domains with sufficient boundary 
regularity. Here, grad^distintr(a^,|/) refers to the Euclidean gradient with respect to x, with 
graddistintr(2;, y) being the gradient with respect to both variables. 

Note also that a consequence of Proposition 12 is that intrinsic geodesies have continuously 
varying directions, and therefore that it makes sense to speak of the angle between a geodesic 
and a Euclidean segment. 

Proposition 14. Suppose that D is a CAT{0) domain, bounded in the Euclidean metric, 
satisfying a uniform exterior sphere condition based on radius r > 0, and a uniform interior 
cone condition based on radius S > and angle a G (0,7r/2]. For every C\ > 0, there exist 
C2, C3 < oo such that, if x,y E D with distintrix, y) > C\ and \/\u — x\ \l \u — x\ < c^, then 

(13) |distj„j,.(u,2/) - (disti„i^(x,y) + |n-x|cos6')| < C3|n-x|^/^ 

where 9 is the angle between the geodesic from y to x and the Euclidean segment from x to 
u that is exterior to the direction from y to x (see Figure 7). Consequently, if x, y (z D with 
X ^ y, then 

(14) grad^ disti„ir(x, y) = -x{.x,y) . 
Moreover, 

(15) graddistintrix^y) = {~xix,y),-xiy,x)). 

Note that Bieske (2010) establishes a similar result for Carnot-Caratheodory spaces. In 
both cases, the relevant distance function satisfies an eikonal equation. 
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u 




Figure 7. Illustration of the configuration of the triangle referred to in the 
statement of Proposition 14. The sides running from y to u and from y to x 
(and of intrinsic lengths dyu and dyx respectively) are intrinsic geodesies. The 
side running from a; to w is a Euclidean segment. 

Proof. In order to demonstrate (13), we establish upper and lower bounds on the difference 

(16) distintr(y, U) - distintr(Z/, x) 

when u is close to x. We abbreviate, setting dy^ = distirAriy^x), etc. 

Let 6' be the exterior angle between the geodesic from x to y and the geodesic from 
X to u. By the CAT(O) property, the Euclidean triangle with the same side lengths as 
a triangle in the intrinsic metric has larger interior angles and therefore smaller exterior 
angles. (The elementary argument for this is given in Bridson and Haefliger, 1999, Chapter 
II. 1, Proposition 1.7(4).) Thus if 6" is the exterior angle of the comparison triangle for x, y 
and u corresponding to the exterior angle 6', then 6" < 6', and so 



dyu = \ dl^ + dl., + 2dyxdxu cos 9" > ^ dl„ + dl,, + 2dyxdxu cos 9' 



yu y ^yx I "'xu ' ^^yx^xu ^^^ " ^ y ^yx I "'xu ' ^^yx^xu 



> \ dl„ + di,, cos^ 9' + 2dyxdxu cos 9' = dyx + dxu cos 9' 



yx I '^xu ' 
> dyx + \u — x\ COS 9' . 

Corollary 13 implies that \9 — 9'\ < c^dxu for small dxu- Hence, 16* — 9'\ < csju — x\ and 
I cos 6* — cos6''| < CsIm — x\. We obtain for |u — x| < C2, for some C2 > 0, 

dyu > dyx + \u — x\ cos 9' > (i,„^. + In — xl cos6' — In — xll cos6' — cos6''| 
(17) y- - y- \ - y- i 

> dyx + I n — x I cos 9 — c^\u — x\ . 
This provides a lower bound on (16) and a bound for one direction of (13). 
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We now establish an upper bound on (16). Fix a point w on the intrinsic geodesic from y 
to X. Then dyx = dy^ + du,x and dyu < dy^ + d^u- We shall require w to be close to x, but 
not as close as u, with |u; — x| = -\/|u — x\ being assumed. 

Because w is close to x and thus also close to u, we may replace the intrinsic geodesies 
from w to u and from w to x by Euclidean segments, without greatly altering lengths and 
segments. Let 9* be the exterior angle at x for the Euclidean triangle with sides d^^w, d^^u 
and du,w From (11), 

(18) \u — w\ < dwu = 1 rP~^l — [i -^ Cq\U — w\ )\u — w\, 

\u — w\ 

(19) \x — w\ < d^x = 1 -\x — w\ < [I + Cq\x — w\ ] \x — w\, 

\x — w\ 

when \u — w\, \x — w\ < 2r sina. 

As before, by Corollary 13, \6 — 6*\ < Cjdxw for small dxw Hence, \9 — 9*\ < c^\w — x\ and 

(20) I cos 6' — cos ^* I < Cii\w — x\. 

These computations allow use to estabhsh an upper bound for distintr(y7 w) ~ distintr(y,a;)- 
First note that 

dyu < dyw + dwu < dyu, + {l + Cq\u — w\'^) \u — w\ 

< dyw + \u — w\+ cq\u — w\^ 

< dyu, + \u — w\+ Cq{\u — x\ + \w — x\Y . 

Now apply the cosine formula to control \u — w|, using (20): 

|u — u;| = y \w — x\'^ + \u — x]"^ + 2\w — x\\u — x\ cos 9* 

{\w — x\ + \u — x\ cos9y + \u — xP sin^ 9 + 2\w — x\\u — x|(cos^* — cos^) 

< \w — x\ + \u — x\cos9 

1 \u — xl"^ sin^ 9 \w — x\\u — x\(cos9* — cos9) 

+ -1 , — , ^ — , — — , — r , 

2\w — x\ + \u — x\cos9 \w — x\ + \u — x\cos9 

1 In — a;Psin^6' lu; — xlln — xl 

< \w — x\ + \u — x\ cos9 H — j j + j j rCslw — x\. 

2\w — x\ — \u — x\ \w — x\ — \u — x\ 

If we take |w — x| = \/\u — x|, with |m — x| < Cg for a suitably small Cg > 0, then 

Im — tfl < Itf — xl + In — xl cos6' + ciolu — xP'^. 
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Combining these bounds implies 

dyu < dy^a + \u - w\ + ce{\u - x\ + \w - x\)^ 

< dyw + \w — x\ + \u — x\ cos9 + cio\u — xf'"^ + cq{\u — x\ + \w — x\)'^ 

< dyw + dw,x + \u — x\ COS 6 + cu\u — xl^'"^ 
= dyx + \u — x\ cos 6 + cu\u — x\'^''^, 

which provides an upper bound on (16). It follows from the above inequality and (17) that 

\dyu — {dyx + \u — x\ cos 6)\ < c^lu — xf''^ , 

which yields the bound in (13). The formula in (14), for the gradient of the intrinsic distance 
distintr(2;, 2/) with respect to x, follows immediately. 

We still need to demonstrate the formula in (15). Let ci, C2 and C3 be as in the statement 
of (13). Fix x,y & D and suppose that distjntr (a^, y) > 2ci. Suppose that u,v G D, 
■\/\u — x| V |m — x| < C2 Aci/4 and \/\v--y\y \v — y\< C2 Aci/4. Let 6x be the exterior angle 
between the geodesic from x to y and the Euclidean segment from x to u. Similarly, let 6y 
be the exterior angle between the geodesic from y to x and the Euclidean segment from y 
to V. Also, let 6' be the exterior angle between the geodesic from y to u and the Euclidean 
segment from y to v. Then by the above reasoning 

(21) |distintr(ii,y) - (distintr(x,y) + |u-a;|cos6'^.)l < csln-a;^/^ 
and 

(22) |distintr(^',M) - (distintr(y,M) + |^^-y|cos6'')| < c^lv - y\^^^ . 

Recall that the Euclidean triangle with the same side lengths as a triangle in intrinsic 
metric has larger interior angles. Using a triangle inequality for angles, \0y — 0'\ is less than 
the angle at the vertex corresponding to y in the Euclidean triangle with sides dxy, dxu and 
dyu- It follows that \0y — 6'\ < Ci2dxu/dxy < Cislu— x| and therefore | cos 6*^— cos 6*' | < Ci3|n— x|. 
This and (22) yield 

(23) \distintr{v, u) - {distintriy, u) + \v-y\ cos 9y)\ < cs\v - y]^/"^ + cis\u - x\\v - y\. 
The triangle inequality applied to the left hand sides of (21) and (23) implies that 

|distintr('i^,M) — (distintr(a^,y) + 1"^ ~ ^1 COS 9y + \u — x\ cos9x)\ 
< cs\u — xl^'"^ + cs\v — yl^'"^ + cislu — x\\v — y\. 

Consequently, graddistintr(a^, y) exists when distintr(a;, y) is viewed as a function of {x,y) G 
{D X D) \ {{u, u) : u & D} and is given by 

graddistintr(a;,y) = {-xix,y),-xiy,x)) . 

D 
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In Proposition 15, we consider solutions of the differential equation dx = xi^iV) dt for 
pursuit and evasion. Proposition 12 established partial regularity for x(x,i/), which does 
not automatically guarantee well-posedness of solutions (as defined below). However, the 
CAT(O) property, together with boundary regularity, will imply well-posedness, even for 
some discontinuous driving paths y. 

Suppose that y{t), t G [0, Ti], is cadlag, of bounded variation on finite intervals, and takes 
values in D. We will say that x(t), t G [0,Ti], is a weak solution to dx = xi.^^y)^'^ if 
x{t) =x(0) + /px(a;(s),2/(s))ds for alHG [0,Ti 



1 • 



Proposition 15. Let D he a CAT{Q) domain satisfying uniform exterior sphere and interior 
cone conditions. For distinct x,y G D, let x{x,y) be the unit tangent vector at x of the 
geodesic from x to y. We consider the differential equation 



(24) 



dx 



X{x,y)dt 



defined in the weak sense for absolutely continuous paths {x{t) : t > 0} in D, driven by 
paths {y{t) : t > 0}, up until the first time that x and y are equal. The problem is well- 
posed, in the sense that solutions x exist, are uniquely determined by initial values x(0), and 
depend continuously on the initial value x(0) and the driving process {y{t) : t > 0} (using 
the uniform distance metric in both cases). 

Proof. The argument is based on the simpler case when the path y is constant in time, 
which we for the moment assume. In this case, existence follows directly from the existence 
of intrinsic geodesies in CAT(O) domains. To show uniqueness, note that, for two solutions 
X and X of (24), since x and x are absolutely continuous and satisfy the differential equation 
weakly, for almost all s, the time-derivatives of x(s) and x(s) must exist and be given by 
x{x{s),y) and x(x(s),y). Exploiting the differentiability of the intrinsic distance given by 
Proposition 14, for x(s) ^ y and x(s) ^ y, one has 



(25) 



-— distintr(a;(s + t),x{s + t)) 
at 



f=0 



-^disti,tr(r(^)(^),rW(^)) 

dt 



i=0 



where F'^^^ F*^^) are unit-speed geodesies running from 



(26) 



— distintr(r(^)(^),rW(^)) 

dt 



x(s), x(s) to y. We will show that 
< 0. 



J t=o 



Consider a Euclidean triangle abc with side lengths satisfying \ab\ = distintr(r^'*''(0), y), \cb\ = 
distintr(f(')(0),|/) and \bc\ = distintr(F(')(0),f ("^(0)). Let z{t) e ah he a point such that 
\z{t) — a\ = t, and let z{t) G 6c be a point such that \z{t) — c| = t. Then Definition 4 implies 
that 



distintr(F(^)(t),F(^)(t)) < \z{t)-I{t)\ < \Z{0)-I 



distintr(F(^)(0),F(^)(0)). 
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This implies (26). It follows that the derivative on the left hand side of (25) is non-positive; 
therefore x = x ii x(0) = x(0), and so uniqueness holds. 

By considering the behaviour over disjoint time intervals, existence and uniqueness follow 
for the case when y is piecewise-constant, in which case the solution curve x is piecewise- 
geodesic. 

We will establish continuous dependence on the initial position x(0) and the driving process 
y, when y is piecewise constant. Suppose that y, y are two piecewise-constant paths in D, 
and X, X solve 

dx = xix,y)dt, dx = xix,y)dt 

for prescribed initial positions a;(0) ^ y{0) and a;(0) ^ y(0). The solutions x, x satisfy the 
differential equations weakly, and therefore, for almost all s, the time-derivatives of x{s) 
and x{s) must exist and are given by x{^{s),y{s)) and x(^{s),y{s)). Arguing as before, for 
x{s) 7^ y{s) and x{s) ^ y{s), we may construct a CAT(O) comparison for the two triangles 
defined by (a) vertices x{s), x{s), y{s) and (b) vertices x{s), y{s), y{s) (see Figure 8). 
Using this comparison, and continuing until either x{t) = y{t) or x{t) = y{t), the function 
distintr(a^, ^) is dominated by its Euclidean counterpart for a two-dimensional quadrilateral 
which is based on a pair of opposing sides of lengths distintr(a^(s),5;(s)) and disti^triyis) ,y{s)) . 
In detail, and using boldface symbols to indicate corresponding Euclidean comparison 
points, we may argue as follows (see Figure 8). Because side- lengths of comparison triangles 
agree, 

distintr(a:(s),y(s)) = |x(s) -y(s)|, 
distintrix {s),y{s)) = |x(s) - y(s) 
distintr(a;(s),x(s)) = |x(s) - x(s) 
distintr(^(s),y(s)) = |x(s) - y(s) 
distintr(y(s),y(s)) = |y(s) - y(s)| . 

Locating x{t) according to distance from x{s) along the intrinsic geodesic from x{s) to y{s), 
and x{t) according to distance from x{s) along the intrinsic geodesic from x{s) to y{s) (and 
locating comparison Euclidean points in the corresponding way), we find 

distintr(a;(s),x(t)) = |x(s) -x(t)| , 

distintr(J(s),x(t)) = |x(s) -x(t)| . 



) 



) 



) 



Now locate the Euclidean point z at the intersection of the Euclidean line segments x(s), y(s) 
and x(t),x(t), and locate z on the intrinsic geodesic from x{s) to y{s) so that 

distintr(J(s), Z) = |x(s) - z| . 
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Figure 8. Illustration of the CAT(O) comparison argument applied to the 
triangles defined by vertices (a) x{s), x{s), y{s) and (b) x{s), y{s), y{s). The 
corresponding Euclidean triangles have vertices marked with boldface symbols. 

Using comparison arguments and the nature of the Euclidean parallelogram x(s)x(s)y(s)y(s) 
we then see that 



distintr(a;(t),x(t)) < distintr(a;(t), 2;) + distintr(^,a;(t)) 

< |x(t) - z| + |z - x(t)| = |x(t) - X(t)| < max {|x(s) - x(s)|, |y(s) - y(s)|} 

= max{distintr(x(s),x(s)),distintr(2/(s),y(s))} . 

This comparison can also be justified by use of Reshetnyak majorization, however we have 
chosen to present an explicit elementary proof. 

We now consider general y in (24). There exist piecewise-constant functions ?/„ converging 
to y uniformly on compact intervals; let x„ be the corresponding solutions to (24), with 
x„(0) = x(0) e D. If \yn{t) - yn,{t)\ < ci for t e [0,T], then \xn{t) - Xm{t)\ < ci by the 
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argument given above. Since the sequence ?/„ is Cauchy in the uniform norm on [0,T], so is 
the sequence x„, which therefore converges to a function x. 

Recall that we are assuming x(0) 7^ 1/(0). Choose fixed ei,e2 > and let T = inf{t > 
: \x{t—) — y{t—)\ < 2£:i}. By part (3) of Proposition 12, there exists (5i > such that, if 
\ui -M2I > ^i, \vi -V2\ > El, \ui -vi\ < 61 and \u2 - V2\ < Si, then \x{ui,U2) -x{'Vi:'V2)\ < 
e2/T. Suppose that n is large enough so that ||/„(t) — y{t)\ < Si A Si for t G [0,T]. Then 
\xn{t) -x{t)\ < Si A El fort e [0,T] and |x(a;(t),y(t)) - x(x„(t),y„(t))| < £2/^. We obtain, 
for t < T, 



< 



x(t)-x(0)- / xixis),y{s))ds 







\x(t) -x„(t)| + 



Xn{t)-x{0)- / x(x„(s),y„(s))ds 

t pt 

x{.Xn{s),yn{s))ds- / x(2;(s),y(s))ds 

^0 

^1 + 0+ / |x(a;„(s),y„(s)) -x(a;(s),|/(s))|ds 



< 

< Si+ 62. 







Since ei, £2 and 5i can be chosen arbitrarily small, we see that x{t) = x(0)+/g x(a;(s), y{s)) ds 
for all t < inf{t > : \x{t—) — y{t—)\ = 0}. Hence, a; is a solution to (24). 

Uniqueness of the solution of (24), for given x(0) and general y, follows by reasoning as in 
(25) and (26). The continuous dependence of solutions on x{0) and y, for general y, follows 
from the above estimates by approximating y by piecewise constant driving processes. D 



3. CAT(O) AND PURSUIT-EVASION 

We consider the Lion and Man problem in a bounded CAT(O) domain D satisfying the 
uniform exterior sphere and interior cone conditions. Alexander et al. (2006) showed that 
e/2-capture, for given e > 0, must occur for the discrete-time variant of this problem. As we 
will see in Section 4, the Lion and Man trajectories x and y will be weak limits of couplings of 
reflected Brownian motions, with drift and small noise, that arise from our capture problem. 

We therefore modify the Alexander et al. (2006) argument to apply to continuous time; 
the modified argument also supplies an explicit upper bound on the capture time. We will 
only need to consider trajectories x and y that are Lipschitz with constant 1. Note that 
Lipschitz trajectories are absolutely continuous, so that the directions dx/dt and dy/dt 
are defined for almost all times t. 
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One can express the trajectories of Lion x and Man y as functions of time t in the following 
differential form: 

dx = xix,y)dt-u^dL'' , 
(27) dy = Hdt-Uydiy. 

Here, H is assumed to be a pre-assigned, time-varying unit length vector generating the 
motion of the Man, x(x, y) generates the motion of the Lion and is defined as in Proposition 
12, for X ^ y, as the unit tangent at x for the corresponding intrinsic geodesic, while 
i^x £ ■A/'x,r and Uy G Afy^r (for r > satisfying the exterior sphere condition of D as given 
in Definition 6) determine the reflection off of the boundary dD. The vector H is assumed 
to be measurable in t; on account of Proposition 12, x is continuous on x ^ y. The terms 
Ur^dL^, respectively UydL^, are differentials arising from Skorokhod transformations and 
are differentials of functions of bounded variation that increase only when x, respectively y, 
belong to dD, and are then directed along an outward-pointing unit normal so as to cancel 
exactly with the outward-pointing component of the drifts x d t, respectively H dt. 

We note that Skorokhod transformations are uniquely defined for a domain satisfying 
uniform exterior sphere and interior cone conditions (Saisho, 1987) (also compare earlier 
results of Lions and Sznitman, 1984), and they then depend continuously on the driving 
processes (using the uniform path metric). In fact, by the definition of x? ^xdL^ vanishes 
identically, while VydV vanishes identically if {H,v) < whenever y G dD. In particular. 
Proposition 15 applies and guarantees the existence of x and its approximation by piecewise- 
geodesic paths for y determined by H. (We include both the Skorokhod transformation 
differentials in (27) as they will both appear in the stochastic version in Section 4.) 

We base our argument on Alexander et al. (2006, Theorem 12). The proof analyzes the 
greedy pursuit strategy arising from the definition of the vector field x? with the Lion always 
directing its motion along the intrinsic geodesic from Lion to Man. The CAT(O) property 
forces the distance between Lion and Man to be non-increasing, and the Man must run 
directly away from the Lion in order to prolong successful evasion. Since the domain is 
bounded, this will, however, not be achievable indefinitely. 

In order to demonstrate the main result in this section. Proposition 17, we will employ 
the following lemma. 

Lemma 16. Under the greedy pursuit strategy described above, in a CAT{0) domain satis- 
fying uniform exterior sphere and interior cone conditions, and at a time t at which Lion 
and Man locations x{t) and y{t) are differentiable in t, 

— dist intr{x(t),y{t)) = -{1 -\y'{t)\ cos a(t)) 

where a{t) is the angle between the Man's velocity y'{t) and the geodesic running from Lion 
to Man. 
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Proof. This follows immediately from the generalization of Gauss' lemma to such domains, 
as was established in Proposition 14. D 

Alternatively, Lemma 16 follows directly from the first variation formula in CAT(O) spaces 
(Bridson and Haefliger (1999), p. 185; Burago et al. (2001), Ex. 4.5.10). 

Proposition 17. Suppose that D is a bounded CAT{0) domain that satisfies a uniform 
exterior sphere condition based on a radius r > and a uniform interior cone condition 
based on a radius 6 > and angle a G (0,7r/2]. Under the greedy pursuit strategy described 
above, there is a positive constant t^ depending only on the diameter of D and e > (and 
not on H in (27)) such that the Lion will come within distance e/2 of the Man before time 
tc, regardless of their starting positions within D. 

Remark. We use e/2 here rather than e, since a further distance e/2 will be required by the 
stochastic part of the argument. 

Proof. This proof follows Alexander et al. (2006), but is modified (a) to account for the 
continuous time context and (b) because we need to derive a specific upper bound tc on the 
time of e/2-capture. Below, we abbreviate by setting £(t) = disti^trix (t),y{t)). 

Let a be the angle defined in Lemma 16. Note that this is defined for almost all times 
t, since the paths x{t), y{t) are Lipschitz and are therefore differentiable for almost all t. 
Evidently, the Lion will have come within e/2 of the Man by time t unless 

(28) /(l-cosa)ds < C{0)-e/2 < diamintr(£') - £/2 . 

Jo 

Now consider the total curvature of the Lion's path. By Proposition 15, the Lion's path is 

uniformly approximated by pursuit paths driven by discretized approximations to the Man's 

path. If x^"'' is the Lion's path driven by a discretized Man's path y^"^', then the Lion's path 

is piecewise-geodesic, with total absolute curvature given by the sum of the exterior angles 

formed at the points that connect the geodesies that occur when x^"-* changes direction. 

CAT(O) comparison bounds then show the total curvature of x^"^ is bounded above by 

^ ' ^disti.tr(xW,y(«))^ ' 

where summation is over the jumps of the discretized path y^"'\ and a^"-* is the exterior angle 
that the jump Ay^"''' contributes to the geodesic running from x^"'^ to y^"''. 

The total curvature of a path is a lower-semicontinuous function of the path (using the 
uniform topology) for CAT(O) spaces. (This is a special case of a CAT(k) result of Karuwan- 
napatana and Maneesawarng, 2007, referred to in Alexander et al., 2010, Theorem 18.) For 
the sake of completeness, we indicate the short proof for the CAT(O) case. Consider a curve 
q of finite length in a CAT(O) space. Its total curvature TC(g) is the supremum of sums of 
exterior angles of piecewise-geodesic curves interpolating q; a CAT(O) comparison argument 
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shows that these sums of exterior angles increase as the interpolating mesh is refined. Let g" 
be a sequence of curves converging uniformly to q. Furthermore, let g"'™ be the piecewise- 
geodesic curve interpolating g" at the points k2~^ for k = 0, 1, . . .. Then, by definition of 
total curvature, 

TC(g"'") / TC(g") as m ^ cx) . 

Bridson and Haefiiger (1999, Chapter II. 3 Proposition 3.3) observe that the CAT(O) prop- 
erty implies that interior angles are continuous functions of their end vertices and upper- 
semicontinuous functions of their centre vertices. This upper-semicontinuity translates into 
lower-semicontinuity for exterior angles, and hence 

limsupTC(g"'™) > TC(g°°''") , 

n— >oo 

where g°°'™ is the uniform limit of g"'™ as n — )■ oo (here we use the CAT(O) property again) 
and is a piecewise-geodesic interpolation of g at the points /c2~™ for k = 0,1,.... Since 
TC(g) = limTC(g°°'™), lower-semicontinuity now follows from 

limsupTC(g") > limsupTC(g"''") > TC(g°^'™) -^ TC(g) as m ^ oo . 

n— >oo n—xx) 

Consequently, the upper bound (29) provides an upper bound on the total absolute 
curvature of the Lion's path in the limit. Bearing in mind the Lipschitz(l) property of 
y, the total absolute curvature r(t) incurred by x between times and t therefore satisfies 

(30) r(t) < r^^d. 

Jo ^s) 

Assume that C{s) > ejl for s <t. By the Cauchy-Schwartz inequality and (28), 

(31) T{t) 
^ Jo t, y Jo 

V^diamintr(-D) - e/2 ■ Vt . 

Next, we follow Alexander et al. (2006) in applying Reshetnyak majorization (Resetnjak, 
1968; see also the telegraphic description in Berestovskij and Nikolaev, 1993, Section 7.4) to 
generate a lower bound on the total absolute curvature of {x{s) : < s < t}. We provide 
details for the sake of completeness. 

We argue as follows. Reshetnyak majorization asserts that for every closed curve ( in D 
(more generally, in any CAT(O) space), one can construct a convex planar set C, bounded by 
a closed unit-speed curve (, and a distance-non-increasing continuous map (p : C -^ D such 
that 0oC = C; moreover, preserves the arc-length distances along (poC^ and C,. Consequently, 
(p restricted to dC will not increase angles and the pre-images under of geodesic segments 
in C, must themselves be Euclidean geodesies (i.e., line segments). 
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By our assumptions about t, the total absolute curvature of {x{s) : < s < t} is finite 
(see (31)). Fix an arbitrarily small 6i G (0,7r/2). It follows from the definitions of length 
and curvature of a path that, for each n, we can approximate the unit-speed curve {x{s) : 
< s < t} by a piecewise-geodesic curve {z{s) : < s < t'} with the following properties: 

- The curve z is parametrized using arc-length. 

- Note that x and y are continuous, so is xi^iH)-, by Proposition 12 (3). Hence, we 
can choose = to < ^i < • • • < ^n = ^ such that the total absolute curvature of 
{x{s) : tj_i < s < ti} is equal to 7r/2 — 5i for all i, with the possible exception of 
i = n. 

- For every z, there exist ti = t^^ <t] < . . . < t"'' = U+i and Si = s'^ < s} < . . . < s™' = 
Sj+i such that z{s-l) = x{tl) and z is geodesic on [sl,sl~^ ], for all i and j. (Notice 
that the curve z is inscribed in the curve x.) 

- The total absolute curvature of {z{s) : Sj_i < s < Sj} is less than 7r/2. In other 
words, the sum (over j) of exterior angles between {z{s) : s]~ < s < s]} and 
{z{s) ■ si < s < s^"*" } at si is less than it/2. (This is a consequence of z being 
inscribed in x and the CAT(O) property.) 

- The difference between the lengths of {z{s) : < s < t'} and {x{s) : < s < t} is 
less than Si. 

Then we have 

(32) total absolute curvature({a;(s) : < s < t}) > f — — (5i j (n — 1) . 

We apply Reshetnyak majorization to the closed curve formed by {z{s) : Sj_i < s < Si} 
and its chord (the geodesic running from z{ti) back to z{ti^i)). Reshetnyak majorization 
guarantees that the total absolute curvature of {z{s) : Sj_i < s < Sj} dominates the curvature 
of its pre-image in the boundary of a convex planar set Cj. Moreover, the perimeter of its 
pre-image in the boundary Cj has length len{{z{s) : Sj_i < s < Si}), while the remainder of 
the boundary of Ci must be a line segment of length distintr(2;(si), z{si-i)). 

The two-dimensional pre-image of {z{s) : Si_i < s < Si} therefore has total curvature 
bound of |. By two-dimensional Euclidean geometry, we can maximize the ratio of the 
length of the pre-image of {z{s) : Sj-i < s < Si} to the length distintr{z{si), z{si-i)) of its 
chord by considering the case of an isoceles right-angled triangle, in which case the ratio is 
-\/2. Accordingly, we obtain the upper bound 

\en{{z{s) : Si_i < s < Sj}) < V2distintT{z{si), z{si^i)) < A/2diamintr(-D) • 

It follows that a portion of the piecewise geodesic curve z which turns no more than | cannot 
have length exceeding v2 times the intrinsic diameter of the region. (Note this is related to 
the Euclidean diameter by Lemma 10.) This implies that we can control the total length of 
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z and thus the total length of x, with 
(33) 
t-bx = len({x(s) : < s < t}) - 5i < len({z(s) : < s < t}) < ^2 diamintr(^) x n . 

Combining Inequalities (32) and (33), we deduce that 
(34) total absolute curvature({a;(s) : < s < t}) > [— — b\\ [n — V) > 

> f--^i)f ^'~^' 1 

- V2 ; VV2diamintr(/^) 

Recall that r(t) = total absolute curvature({x(s) : < s < t}) and len({x(s) : < s < 
t}) =t. Letting 5i -^ in (34), it follows that 

r(t) TT f 1 1 



t 2 Vv^diamintr(^) ^ 

In combination with (31), this yields 



n 



i v2,,l...M) -l) ^ ^Vd.am,.(D)-./2.i= 



2 V v2diamintr(-D) t 
and hence the quadratic inequality for q = y/t, 

V2v^dianw(I?)y V e ^ '°*^^ ^ ' J "^ 2 - 

The left hand side is negative for q = and the coefficient of g^ is positive, so there is 
exactly one positive root qc (which can be written out explicitly in terms of diamintr(-D) and 
e). Combining this with our earlier arguments, it follows that the Lion will come within e/2 
of the Man by time tc '■= q"^- D 

4. From Brownian shy couplings to deterministic pursuit problems 

This section is devoted to the proof of Theorem 1. Consider a co-adapted coupling of 
reflecting Brownian motions X and Y in the bounded domain D CW^ satisfying uniform ex- 
terior sphere and interior cone conditions. Saisho (1987) showed that the reflected Brownian 
motions can be realized by means of a Skorokhod transformation as strong solutions of 
stochastic differential equations driven by free Brownian motions. As discussed in Section 
1.2, we can use arguments embedded in the folklore of stochastic calculus, and stated expli- 
citly in Emery (2005) and in Kendall (2009, Lemma 6), to represent this coupling as 

(35) dX = dB-uxdL^, 

(36) dY = (j^dE + K^dA) -z/ydL^, 
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where A and B are independent d-dimensional Brownian motions, and JJ, K are predictable 
{d X (i)-niatrix processes such that 

(37) J^ J + K^K = (dxd) identity matrix . 

Here L^ and L^ are the local times of X and Y on the boundary. 

The advantage of this explicit representation of the coupling is that we can track what 
happens to X and Y when we modify the Brownian motion B by adding a drift. We will see 
that the effect of adding a very heavy drift based on the vectorfield xi^y ^) will be to convert 
(35)-(36) into a stochastic approximation of the deterministic Lion and Man pursuit-evasion 
equations (27) over a short time-scale. 

Proposition 18. Suppose that D cM."^ is CAT{0), is bounded in the Euclidean metric, and 
satisfies a uniform exterior sphere condition and uniform interior cone condition. For any 
e > and X and Y satisfying (35)-(36) with ^(0), ^(0) G D, there exists t > such that 



(38) P 



sup distintr{X{s),Y{s)) <e 

t/2<s<t 



> 0. 



Proof. Consider the following modification of (35)- (36), 
(39) 

X"(t) = X{0) + B{t)+ [ nx{X^{s),Y^{s))ds- [ z/x"(.)dLf , 

Jo Jo 

(40) 

y"(t) = F(0)+ f (jJd5(s) + KjdA(s))+ / njJx(^"(s),F"(s))ds- / i/y.(,)dL 
Jo Jo Jo 

By the Cameron-Martin-Girsanov theorem, the distributions of the solutions of (35)-(36) 

and (39)-(40) are mutually absolutely continuous on every fixed finite interval. We will show 

below that, after rescaling time, paths of {X"'{-),Y"'{-)), for large n, will be uniformly close 

to those for the corresponding Lion and Man problem. Application of Proposition 17 will 

then enable us to finish the proof. 

We will make the following substitutions, X"(t) = X"(nt), y"(t) = r"(nt), B{t) = 

B''int)/y^, A{t) = A''{nt)/y^, J(t) = |(")(nt), K{t) = KW(nt). Then (39)-(40) take the 

form 

(41) x^(t)=X{0) + ^B-{t)+ [ x{X-{s),Y^{s))ds- [ ^^„,,,dLf , 



yn 



n .In ./Q 



(42) r"(t) = r(o) + ^ [ f(j("yd5"(s) + (K("ydA"(s 

+ [\lryx{X^{s),Y^{s))ds- f Uy^,dL^^ 
Jo Jo 



SHY COUPLINGS, CAT(O) SPACES, AND THE LION AND MAN 32 

Note that i?" and A^ are Brownian motions. 

Consider the analog of (41)-(42), but without boundary: 

(43) ^"(t) = ^5"(t) + [ x(^"(s),F"(s))ds, 
'^ Jo 



^ 'JO 



(44) \/"(t) = ^ / ((jW)Td5'^(s) + (K("))^dA"(s)) + / (J("))^x(^"(s),l""(s))ds 



All components of the sextuplet 

(45) 

W"(t) = (^"(t),^5"(t), f x(^"(s),?"(s))ds, 

V-{t),^j ((J("))^d5"(s) + (KW)TdI"(s)),y"(Ji"))Tx(^"(s),?"(s))ds 

are tight by the criterion given by Stroock and Varadhan (1979, §1.4) since the diffusion coef- 
ficients and the drifts are bounded by 1. So, on an appropriate subsequence, W" converges 
weakly to a limiting process W°°. By abuse of notation, we will denote this subsequence 
W". In particular, U^{t) and V^lt) converge weakly, so, by Saisho (1987, Thm. 4.1) (which 
applies because of the conditions imposed on D), (X",y") converges weakly to a limiting 
continuous process {X°°,Y°°) along the same subsequence. It follows that 

(46) 

Z"(t)= (X"(t),y"(t),t/"(t),^5"(t), fxiX^{s),Y^is))ds, 

\/"(t),^ f (^(|("))Td5"(s) + (K("y dA"(s)) , f (j("))Tx(^"(s),r"(s))ds 



is tight and, therefore, converges weakly along a subsequence. Once again, we will abuse the 
notation and assume that Z" converges weakly. By the Skorokhod Lemma (Ethier and Kurtz, 
1986, Sec. 3.1, Thm. 1.8) we can assume that the sequence Z" converges a.s., uniformly on 
compact intervals. 

The fourth and seventh components of Z" are Brownian motions run at rate - so they 
converge to the zero process as ra — )■ oo. The fifth and eighth components of Z"' are both 
Lip(l); their limits are therefore also Lip(l). These observations and (43)-(44) imply that 
the limits V°° and U°° of V"" and f/" are Lip(l). 
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Let T* = mf{^> : X^{t) = F°^(t)}. The bounded vectorfield x(^",>'") depends 
continuously on X" and F" (Proposition 12). We may therefore apply the dominated con- 
vergence theorem and (43) to deduce the following integral representation for U°° , 

(47) f7°°(t)= /" x(^°°(s),?°°(s))ds fort<f*. 



Recall that, by the Skorokhod representation, we can assume that X"(t) and y"(t) con- 
verge almost surely. Lemma 19 proved below shows that X°°, Y°° are both still Lip(l), with 
respect to the intrinsic metric of D. Hence we can apply the results on CAT(O) Lion and 
Man problems at the end of Section 2. 

Fix an arbitrarily small e > 0. It follows from (47) and from Proposition 17 that there 
exists ti < oo not depending on X(0),F(0) or cj, such that distintr(^°^(i), ?°°(i)) < £^/2 for 
t > ti. We conclude that for some no < oo, depending on X{0) and Y[0), and all n > no, 



P 



sup distintr(^"(t),F"(t))<£ 
ti<t<2ti 



>0. 



Changing the clock to the original pace, we obtain 



P 



sup distintr(X"(t),F^Xi))<^ 

ti/n<t<2ti/n 



>0. 



By the Cameron-Martin-Girsanov theorem. 



(48) 



P 



sup dist,,,triX{t),Y{t))<e 

ti/n<t<2ti/n 



>0. 



Since e > is arbitrary, this completes the proof. 



D 



Lemma 19. Let D be a domain satisfying uniform exterior sphere and interior cone condi- 
tions. Suppose that Z is a continuous process on D derived by the Skorokhod transformation 
from a free process S that has Lip{l) sample paths. Then Z itself has Lip{l) sample paths 
with respect to the intrinsic metric. 

Proof. Following Saisho (1987, §3), consider the step function Sm obtained from S by 
sampling at instants k2~"^, for A; = 0, 1, .... Suppose that 2"™ < r, where r is the ra- 
dius on which the uniform exterior sphere condition is based. Let 1 be the projection onto 
D described in Lemma 7. The Skorokhod transformation of 5*^ is Z^, given by projecting 
increments back onto D: 



(49) Z^t) 



Zmiik - 1)2—) + ASm{k2-"^) 
Z(0) 



for k2-"' <t <{k + 1)2' 
for < t < 2-™ . 



On account of the Lip(l) property of S, this projection is defined when 2 



< r. 
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From Saisho (1987, Theorem 4.1), we know that Z^ — ^ "Z uniformly on bounded time 
intervals. We compute the maximum possible Euclidean distance between Z^is) and Zm(t), 
if < t — s < 2~™, when one or both of 2'"s, 2"^t are non-negative integers. Since Z^ is 
constant on intervals [k2~"^, [k + 1)2~™), it suffices to produce an argument for the case 
when 2™s = k — 1 and 2"^t = k. We therefore proceed to bound the Euclidean distance 
\Zm{k2^'^) — Zm{{k — 1)2"™) I . We will show that this can only exceed 2^™ by an amount 
which, for large m, will make a negligible contribution to path length when summed over 
the whole path. 

If Zm{k2~'"^) ^ dD, then there is nothing to prove, since the jump is ASm{k2~''^), which 
is bounded in length by 2^"* since S is Lip(l). So we instead suppose that Zm{k2~"^) G dD. 
For convenience, set y = Z^,{k2-"') - {Z^{{k - 1)2""^) + ASm{k2-"')) to be the Skorokhod 
correction to be applied at this step, and set a = \ASm{k2~"^)\ to be the length of the 
uncorrected jump. Finally, let 6 be the angle between the vector y and the negative jump 
—ASm{k2~'^). These definitions are illustrated in Figure 9, together with the supporting 
ball B at Zm{k2~"^) G dD whose centre is located at Zm{k2~'^) — Xy for some X = r/\y\ > 
and whose existence is guaranteed by the construction of the x i— )■ x projection map as 
described in Lemma 7. 




ZJk2- 



ZJ(k-l)2-^) 



"m 



Figure 9. Illustration of the geometry underlying the argument of Lemma 19. 

First note that \Zm{k2~'^) — Zm{{k — 1)2~"^)| = \/a? + y"^ — 2ay cos 6 (where we abuse 
notation by letting y also stand for the length of the vector y). This increases as 6 increases 
to 71, so long as a, y, Zm{k2~"^) are held fixed. Thus we can assume that 9 has increased to 
the point where Zm{{k — 1)2^™), as well as Zm(A;2~™'), belong to dB. (This will happen if, 
as required above, 2""^ < r.) Now observe that the distance \Zm{k2^"^') — Zm{{k — 1)2~™)| 
will be bounded above by the smaller of the two distances from Zm{k2~"^) to the intercepts 
of dB by a line parallel to y, and at distance a from Zm{k2~'^). But two applications of 
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Pythagoras' theorem show that this distance is given by 



35 



^J a? + (r — ^/r'^ — a?Y 




for some z* in the range [0,^]. (The last step arises from a second-order Taylor series 
expansion.) Therefore 



y a? + (r — ^/r"^ — aP'Y 



< ^'•■/i^O + i^) 



< 



-ii + ^ 



(using a/1 + z < 1 + 1^ for ^ > —1). 
Thus the total path length over the time interval (s, t) is bounded above by 



which converges to t 



((t-s)2'" + 2) x2- 



s as ?7i — 7- oo. Hence we obtain 

distintr{Z {s),Z{t)) < t-s 



D 



thus establishing the Lip(l) property in intrinsic metric for Z. 

We will show that the bound in Proposition 18 is uniform over all X(0) and 1^(0). We 
will switch from the intrinsic distance to the Euclidean distance in the formulation of the 
next proposition. This is legitimate in view of (10). 

Proposition 20. Let D be a domain satisfying uniform exterior sphere and interior cone 
conditions. Suppose that there exist ti > and £i > such that, for any X and Y satisfying 
(35)-(36) with X{0), Y{0) e D, 



P 



inf \X(t)-Y(t)\ <ei 

0<t<ti 



> 







inf \X(t)-Y(t)\ <£i 

0<t<ti 



> 



Pi 



(50) 

Then 

(51) P 

for some pi > not depending on X{0) and Y{0). 

Proof. Suppose (51) does not hold. Then there exist ti > 0, £i > 0, sequences {x„}„>i, 
{yn}n>i of points in D, random processes {At,t > 0}, {Bt,t > 0}, {JI",t > 0} and {K^,t > 
0}, and solutions X" and Y"' of (35)-(36) satisfying the following properties. The processes 
A and B are (i-dimensional Brownian motions starting from 0, and independent of each 
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other. The {d x (i)-matrix-valued processes J" and K" are predictable with respect to the 
natural filtration of A and B, such that (J")''' J" + (K")'''Kt is the {d x d) identity matrix at 
all times t. Let X"- and Y"" denote solutions to (35)- (36) based on the Brownian motions A 
and B, using the predictable integrators J" and K*^, and starting from X"(0) = Xn ^ D and 
r"(0) =y^eD. Then 



(52) P 



inf |X"(t)-F"(t)| >ei 
<t<ti 



> 



Let {Mt'^.Mt''^) = (/Jd5„/p(J^)^d5, + J^iJ^ydA,). The processes iVf^'^ and M""'^ 
are Brownian motions and so the sequence of pairs is tight, which therefore possesses a 
subsequence converging in distribution. By abuse of notation, we assume that the whole 
sequence (M"'\M"'^) converges in distribution to, say, (M°°'\ M°°'^). It is clear that M°°'^ 
and M°°'^ are Brownian motions. 

Let I't = cr((Mf''\ Mf='2), s < t) be the natural filtration for (M°°'\ M°°'2). We will show 
that (M°°'^,M°°'^) are co-adapted Brownian motions relative to {J^t}- Since (M"''^,M"'^) 
are co-adapted Brownian motions, for all < ti < ^2 £ ■ ■ ■ ^ ^n ^ ^ ^ si < S2, the random 
variable M"^'^ — M"^'^ is independent of 

{{M-'\ M,f ), (M^/, M,f ), . . . , {M^\ M^f)) . 

Independence is preserved by weak limits, so M^'^ — M^'^ is independent of 

((M,7'\ M,7'^), (M,^'\ M,^'^), . . . , {Mrj\ Mr/)) . 

This implies that M^'^ — M^'^ is independent of J-^. Since the same argument applies to 
M~'2-M,^'2, we see that (M°°'\ M^'^) are co-adapted relative to {Ft}. Recall from Section 
1.2 that this implies that there exist Brownian motions {Af^,t > 0} and {5j°°,t > 0} and 
processes {J^, t > 0} and {K^,t > 0} such that iM^'\ M^'^) = (/J d 5f , J^OTV d 5,°° + 

/J(Kr)^dAr). 

Recall that (M"'\ M"'^) — )■ (M°°'\ M°°'^) weakly in the uniform topology on all compact 
intervals. Going back to the original notation, we see that 

(J dB,, J iJ^:)^ dBs+ f (K^y dA?\ ^ 

di?r, / (JD'^dEr + / i^TVdAT 

Jo Jo 

weakly in the uniform topology on all compact intervals. By the Skorokhod Lemma, we can 
assume that the processes converge a.s. in the supremum topology on compact intervals. 

Since D is compact, we can assume, passing to a subsequence if necessary, that the initial 
points satisfy x„ — )■ Xqo € D and y„ — )■ yoo G -D as n — )■ cxd. In view of the representa- 
tion of coupled reflected Brownian motions using stochastic differential equations (35)-(36), 
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established in Saisho (1987, Theorem 4.1), and employing the continuous dependence on 
driving Brownian motions established there, we see that (X", y") — )■ {X°°,Y°°) weakly in 
the uniform topology on all compact intervals, where {X°°,Y°°) represents the solution to 
(35)-(36) with X°°(0) = x^o, F°°(0) = yoo, corresponding to A°^, 5~, J°° and K°°. We 
obtain from (52) and weak convergence of {X"-,Y"-) to {X°°,Y°°) that, for every n, 



P 



inf |X°°(t)-r°°(t)| >ei 

0<t<ti 



> 



Taking the limit as ra — )■ oo, this contradicts (50) in the statement of the Proposition. 
Consequently (51) must hold for some pi. D 

We now complete the proof of Theorem 1, applying Proposition 20 together with stand- 
ard reasoning. Consider processes X and Y starting from any pair of points in D and 
corresponding to any "strategy" J and K. Because of the uniform bound in Proposition 
20, the probability of X and Y not coming within distance Si of each other on the interval 
[kti, {k + l)ti], conditional on not coming within this distance before kti, is bounded above 
by 1 —pi for any k, by the Markov property. Hence, the probability of X and Y not coming 
within distance ei of each other on the interval [0, kti] is bounded above by (l—pi)^. Letting 
k -^ oo, it follows that X and Y are not ei-shy. Since ei can be taken arbitrarily small, the 
proof of Theorem 1 is complete. 

We remark that the matrices J and IK employed in (35)- (36) are predictable and, con- 
sequently, the choice of the pursuer's velocity is made without looking into the future. This 
is in contrast to the pursuit-evasion problems and associated paradoxes discussed by Bollo- 
bas. Leader, and Walters (2011). 

5. Complements and conclusions 

We conclude this paper by remarking on some supplementary results and concepts, and 
by considering possibilities for future work. 

5.1. Comparison with previous methods. The fundamental idea in this paper turns 
out in the end to resemble that of Benjamini et al. (2007), but uses simple notions of weak 
convergence and tightness, rather than detailed large deviation estimates. Moreover, the use 
of metric geometry notions enables us to finesse many analytical technicalities. (Perhaps 
this is the first application of modern metric geometry to Euclidean stochastic calculus?) 
On the other hand, the stochastic control methods of Kendall (2009) are quite different. 
The stochastic control approach uses potential theory to estimate the value function of an 
associated stochastic game; consequently the methods of Kendall (2009) may be expected to 
give sharper information (bounds on expectation of stopping times), but in more limited cases 
(convexity of domain). However, one can observe that, at least in principle, the stochastic 
game formulation still applies in the general case. For example, there is a value function to 
be discovered for a stochastic control reformulation of Theorem 1, and in principle it might 
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be possible to estimate this value function and so gain more information than is supplied by 
the weak geometric bounds established above. 

We note that many promising ideas based on stochastic calculus fail to show non-shyness 
because they cannot be applied to "perverse" couplings with the property that, on some 
time intervals, |X — y| grows at a deterministic rate (see Example 4.2 of Benjamini et al., 
2007). 

Also note that the proof in Kendall (2009), which works in convex domains, does not 
appear to be (directly) extendable to calculations involving the intrinsic metric - simple 
manipulation using symbolic Ito calculus (Kendall, 2001) shows that the drift of distintr(-^, Y) 
is unbounded at distances bounded away from zero. In particular, Bessel-like divergences for 
distintr(-^, y) of magnitude a occur when the geodesic from X to y touches a concave part 
of dD at X and |x — y | = 1/a. The first-order differential geometry given in Proposition 14 
(the generalization of Gauss' lemma) is the best we can do for CAT(O) domains satisfying 
uniform exterior sphere and interior cone conditions. 

5.2. Higher dimensions and the failure of CAT(O). For planar domains, CAT(O) and 
simple-connectedness are equivalent, in which case, by Theorem 2, there are no shy co- 
adapted couplings. In higher dimensions, it is natural to ask whether the CAT(O) condition 
is essential for there to be no shy coupling. We do not at all beheve this to be the case. It 
is possible to give an argument suggesting that star-shaped domains with smooth boundary 
conditions cannot support shy couplings, by establishing the analogous result for a corres- 
ponding deterministic pursuit-evasion problem. To apply this argument to the probabilistic 
case would require more careful arguments. We therefore leave this as a project for another 
day. 

As a spur to future work, we formulate a bold and possibly rash conjecture: 

Conjecture 1. There can be no shy co-adapted coupling for reflecting Brownian motions in 
bounded contractible domains in any dimension. 

While resolution of the star-shaped case appears to be largely a technical matter, we 
believe that new ideas will be required to make substantial progress towards resolving the 
conjecture. 

5.3. When can shyness exist? Many examples of shy couplings can be generated using 
suitable symmetries. However, we do not know of any examples in which symmetries play 
no role. Accordingly we formulate a further conjecture: 

Conjecture 2. // a bounded domain D supports a shy co-adapted coupling for reflecting 
Brownian motions, then there exists a shy co-adapted coupling that can be realized using a 
rigid-motion symmetry of the domain D. 

A stronger from of the above conjecture, saying "If a bounded domain D supports a shy 
co-adapted coupling for reflecting Brownian motions, then the shy coupling is realized using 
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a rigid-motion symmetry of the domain D," is false. To see this, consider the planar annulus 
A = B{0, 2) \ i3(0, 1) and let T be the symmetry with respect the origin. Let X be reflected 
Brownian motion in A and Y = T{X). Let D = A x (0, 1) and let Z be reflected Brownian 
motion in (0, 1), independent of X and Y . Then (X, Z) and (F, Z) form a shy coupling in 
D which cannot be realized using a rigid-motion symmetry of D. 

Note that Benjamini et al. (2007, Example 3.9) supplies an example based on Brownian 
motion on graphs, for which there is no fixed-point-free isometry and yet a shy coupling 
exists. However we do not see how to use the idea of this construction to construct a 
counterexample to the above conjecture. 

5.4. Further questions. We enumerate a short list of additional questions. 

(1) Shyness is interesting for foundational reasons: coupling is an important tool in prob- 
ability, and shyness informs us about coupling. We do not know of any honest applica- 
tions of shyness. However, one can contrive a kind of cryptographic context. Suppose 
one wishes to mimic a target Y , which is a randomly evolving high-dimensional struc- 
ture, in such a way that the mimic X never comes within a certain distance of the 
target Y . Shyness concerns the question, whether it is possible to do this in a way 
that is perfectly concealed from an observer watching the mimic X alone. 

(2) In this formulation, it is not clear why one should restrict consideration to co-adapted 
couplings. Our methods do not lend themselves to the non-co-adapted case, and the 
question is open whether or not results change substantially if one is allowed to use 
such couplings. In particular, it seems possible that Conjecture 2 might have a quite 
different answer in this context. 

(3) In further work (Bramson et al., 2011) we plan to study the deterministic pursuit- 
evasion problem, in conjunction with shy couplings, for multidimensional CAT(k) 
domains possessing "stable rubber bands," a condition that is partly topological and 
partly geometric. As a corollary, we plan to prove that there are no shy couplings in 
multidimensional star-shaped domains. 

(4) The Lion and Man problem has been generalized to the case of multiple Lions. (An 
early instance is given in Croft, 1964.) Can one formulate and prove useful results 
for a corresponding notion of multiple shyness? 
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